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Abstract 

In the framework of the seesaw models with triplets of fermions, we evaluate 
the decay rates of \x — > and r — > Z7 transitions. We show that although, due to 
neutrino mass constraints, those rates are in general expected to be well under the 
present experimental limits, this is not necessarily always the case. Interestingly 
enough, the observation of one of those decays in planned experiments would 
nevertheless contradict bounds stemming from present experimental limits on 
the \x — > eee and r — > 31 decay rates, as well as from \i to e conversion in atomic 
nuclei. Such detection of radiative decays would therefore imply that there exist 
sources of lepton flavour violation not associated to triplet fermions. 
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1 Introduction 



The search for flavour changing rare leptonic decays, in particular for [i — > cy, r — > fi'j 
and r — > cy decays, has been the object of intense experimental investigations for 
decades [I]. With respect to the present experimental upper limit, Br(/i — * ej) < 
1.2- lO^ 11 11], Br(r ^7) < 4.5 -lO" 8 [2], Br(r cy) < 1.1 • 10^ 7 [TJ, new experiments 
are expected to improve in the near future their branching ratios by as much as three 
orders of magnitudes for the first decay mode [3] and by one or two for the two others [I] . 

The recent experimental evidence for neutrino masses has shown that lepton flavour 
is violated in the neutrino sector and that, consequently, in a model independent way, 
these decay rates are predicted to be different from zero. The actual predicted rate, 
however, turns out to be highly model dependent. There are three basic models which 
can explain the neutrino masses at tree level, from the exchange of heavy states, through 
the seesaw mechanism. The above rare decays have been studied at length in the frame- 
work of two of these seesaw models, with right-handed neutrinos [5] (type-I seesaw [Fj]) 
and with one or several Higgs triplets [7j (type-II seesaw [8]). In this letter we perform 
the calculation of these decay rates in the framework of the third seesaw model, with 
heavy triplets of fermions (type-Ill seesaw j9]). This model has been studied in detail, 
both from the theoretical and phenomenological point of view, in Ref. [TU], where the 
result on these rare decays has already been presented without the detailed calculation. 
This letter also contains a determination of the constraint that \i to e conversion in 
atomic nuclei implies on the type-Ill seesaw model. 



2 The type-Ill seesaw Lagrangian 

The type-Ill seesaw model consists in the addition to the standard model of SU(2) 
triplets of fermions with zero hypercharge, S. In this model at least two such triplets 
are necessary in order to have two non-vanishing neutrino masses. A non-vanishing 
h hi rate can nevertheless be induced already with only one fermionic triplet. In 
the following, we will not specify the number of triplets so that our calculation is valid 
for any number of them. Being in the adjoint representation of the electroweak group, 
the Majorana mass term of such triplets is gauge invariant. In terms of the usual and 
compact two-by-two notation for triplets, the beyond the SM interactions are described 
by the Lagrangian (with implicit flavour summation): 

£ = Tr[Sip£] - -Tr^M^ + I>M*£] - ^Ey/TY^L - Lv^l^E^ , (1) 
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with L = (I, v) T , <p = (0+ 0°) T = (0+ (v + H + ir])/V2) T , 4> = iT 2 </>*, S c = Cjf and 
with, for each fermionic triplet, 

TP/y/2 S+ \ ( £ 0c /V2 S~ c 

E- -S° /y/2 J ' V S+c -E 0c /V2 

Without loss of generality, in the following we will assume that we start from the basis 
where is real and diagonal. In order to consider the mixing of the triplets with the 
charged leptons, it is convenient to express the four degrees of freedom of each charged 
triplet in terms of a single Dirac spinor: 

^£+ c + £ fl . (3) 

The neutral fermionic triplet components on the other hand can be left in two-component 
notation, since they have only two degrees of freedom and mix with neutrinos, which 
are also described by two-component fields. This leads to the Lagrangian 

C = ^i0> + Y^ipY? R - *M E * - fe^££ + h.c. 

+ g {w+W Rl ,P R m + W+V% 7m P l ¥ + h.c.) - g 

(fW R Y^ L + V5/W S I L + <P + W R Yd L - y/2<j> + u?Y£V + h.c.) . (4) 
The mass term of the charged sector shows then the usual aspect for Dirac particles: 

£ H«(s:)(^)-ra(lt)(^)^ < 5 > 

with v = y/2((j) ) = 246 GeV. The symmetric mass matrix for the neutral states is on 
the other hand given by 



C 3 -(FF£ 0c 



Y^v/2V2 \ 
Y s *v/2V2 M E /2 J 



vo 



Y s v/2V2 M E /2 ) { £ 0c 



(6) 



Diagonalization of the mass matrices 

To calculate the l\ — > I2J decay rates, we will work in the mass eigenstates basis. As 
it happens with any Dirac mass, the charged lepton mass matrix can be diagonalized 
by a bi-unitary transformation 

'&)-«*•(&)• (7) 
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where Ul,r are (3+n)-by-(3+n) matrices, if n triplets are present. On the contrary, the 
symmetric neutral lepton mass matrix can be diagonalized by a single unitary matrix 



Vh \ - 11 ( V ' L 

s o c J - Uo y S /0c 

It is convenient to write the mixing matrices in terms of three-leptons-plus-n-triplets 
sub-blocks 



In the following we will calculate the decay rates at 0((Yj^v j My) 2 ) , which is a good 
approximation as long as M% is sufficiently big compared to Y^v. In order to do so it 
can be checked that it is enough to calculate all the mixing matrix elements at order 
0([(Yxv,mi)/Mx} 2 ). We obtain: 



e' 



U Lll = l-e U LW = Y^M^v U Lm ^-M^Y s v U L y* = 1 

U Rll = 1 Urm = miY^M^v U Rm = -M^ 2 Y s miv U R ^ = 1 

Uq uu = (1 — ^)UpMNS U Q vY; = ^E^E 1 ^ Uqy,u = -^s^E^^Ow f^OEE = (1 — \) 

(10) 

where e = yY"£M s 2 Y" s , e' = y^s^s^s^s 1 an< ^ Upmns is the lowest order neutrino 
mixing matrix which is unitary. Note that e is nothing but the coefficient of the 
unique low energy dimension-six operator induced by the triplets, once they have been 
integrated out [10] Eq. ffTUj) shows as expected that the (3 + n)-by-(3 + n) mixing 
matrices Ul,r.o are unitary but the various submatrices are not. The neutrino mass 
matrix in this model is given b}0: 

V 2 rr 1 , , 

m u = Yl Ye- (H) 

2 E My, K J 



Lagrangian in the mass basis 

After the diagonalization of the mass matrices, we obtain the following Lagrangian 
in the mass basis (omitting from now on the primes on the mass eigenstate fields): 

£ = C-Kin + £-CC + C-NC + C-H,r) + £<[>- , (12) 

lr rhe e' contribution does not appear in the low energy effective theory as it involves external S's. 

2 As for the masses of the charged leptons, they are essentially unaffected by the presence of the 
E's as the difference between the physical masses of the I' and the ones of the Z's, mi, is of order 
miYlv 2 /Ml. 
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where 



(13) 



C CC = j= ( 7 * ) y (p L g?> + Prf^) ( £ ) + h.c. 

c *° - ( 1 * > 7 " z " (/5i9 " c + prfC) ( J ) (14) 

+ ! 2^( ? *»"<^ + F « 9 «»(*) <15) 



with 



q cc _ ( 9% 9% \ _ (l + e)U „„ -Y^f 2 
~ U?£ 9 C L l)-{ ^2(1-1) 

9 C rZ 9<k% \_ ( - mi YlM^v 



(17) 



9L — „NC „NC I — 1 «f-lv „, ™„2/] l ly J 



" 1 9& 9 N nt ) ~ V M" 2 W -cosH w ) (20) 

flg, 9L\-( m,(3e-l) -m^M^ \ 

~ 1 9 H Rm ) ~ { -M^m lV ... ) (22) 

9l u 9l l9 \_ ( - (e + 1) m, (1 - e) - mj^M^ \ 

" 1 7 " V -YW (1 - e) + M E - 2 F s m 2 , ... J ^ 



and 



9l=mU^ v { 9t=mYlM^^ 2 

9% = - (1 - e) \ 4 = (1 - 6) ^ (1 - ^ ) - v^^M^ 



.(25) 



The dots in Eqs. (l2~Tj) - ([24"j) refer to \E'-\I/ interactions which we omit here since they do 
not contribute to the one-loop l\ — > l-fi rates. 
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Figure 1: Diagrams contributing to [i — > e7- , 77 are the three Goldstone boson 
associated with the and Z bosons. H stands for the physical Higgs boson. 

3 fi — ► e7 and r — ► Z7 decays 

In the following we perform the calculation of the \x — > &-/ rate. The r decay rates 
will be obtained straightforwardly from it later on. As it is well-known, the on-shell 
transition fi —> is a magnetic transition so that its amplitude can be written, in the 
m e — > limit, as : 

T - erf) = A x W e (p - q) [iq u e x a Xu (1 + Ts )] u M (p) , (26) 

with e the polarization of the photon, p^ the momentum of the incoming muon, the 
momentum of the outgoing photon and a^ v = | [7^,7*,]. Using the Gordon decompo- 
sition we can rewrite it as 

T(^e7) = 4x«;(p- 9 )(1 + 75) (2p-e- m^) (p) . (27) 

In the following we will calculate only the p ■ e terms. The terms proportional to ^ can 
be recovered from the p ■ e terms through Eq. ( 1271) . All in all, this gives: 

•* 

Tft 

r(/i-e 7 ) = _qA| 2 . (28) 

fi — ► e7 amplitude and decay rate 

In the mass eigenstate basis, from the Lagrangian of Eqs. (II 3ft - (TIB"]) , there are four- 
teen diagrams contributing to [/, — > ej, as shown in Fig. 1. The detailed calculation is 



5 



presented in the appendix | . In the limit in which M s > M w , at 0((^) 2 ), the total 
amplitude is given by: 



.Gp M e 



^m^Ue {p - q) (1 + 75) i(T X uS X q v u^ (p) 



X 



t 

PAINS 



(29) 



where C = —6, 56 and % v . = -rrt- Note that the second term is the usual contribution 
from neutrino mixing [12], while the first one is the explicit contribution of the fermion 
triplet (s). As well known, a GIM cancellation operates in the second term. The total 
decay rate is then given by: 



r (ji -> erf) 



G s F M2 e 2 ml 
8192tt 5 

and the branching ratio reads 



PMNS) e i [V PAINS 



'I 1 



Br (fi -> e 7 ) = — ^ 



+ Cj e eyU - y^z^ (UpMNs) ei (Upmns 



tft 



(30) 



(31) 



r — > Z7 decays can be obtained from Eq. (1311) by replacing /x by r, e by / and by 
multiplying the obtained result by _Br(r — > eu T u e ) = (17.84 ± 0.05) • 10~ 2 [T]. 



4 Phenomenology 

From the result above, it is not surprising that in general we expect a very tiny \x — > 
rate. For instance, omitting flavour indices, for a given value of Ms we would ex- 
pect in general from the seesaw formula that Y£ ~ m u Mz/v 2 ~ \J 5rn 2 atm MY< / v 2 ~ 
M s /(10 15 GeV). This gives e ~ m„/M s ~ 10~ 25 (10 15 GeV/M E ) and x v ~ 5m 2 atm /M^ v ~ 
10" 24 which leads to Br (fi — > ej) ~ 10~ 52 • (10 15 GeV/M s ) 2 , far below the present up- 
per limit 1.2- 10 -11 . In this case, even for as low as 100 GeV, we get Br (/i — > ej) ~ 
10 -26 . Similarly, for r — > /X7 and r — > &y, we get both rates of order 10 _53 (10 15 GeV/Ms) 2 , 
far below the present upper limit 4.5 • 10~ 8 and 1.1 • 10~ 7 , respectively. 

There are cases, however, in which the branching ratio can be much larger without 
any fine-tuning of the Yukawa couplings and mass parameters. This is the case if 
neutrino masses are generated through "direct lepton violation" (DLV) (see Ref. [TU]). 

3 General formulae for radiative fermion decays have been derived in detail in Ref. [IT], although 
restricted to the case in which all fermion masses arise from the standard Higgs mechanism. In 
consequence, isospin invariant mass terms as those essential in seesaw models were not taken into 
account. 
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i.e. if neutrino masses are directly proportional to a small lepton number violating 
scale rather than inversely proportional to a high scale. DLV appears naturally in the 
type-II seesaw model, since two scales are present there: the mass of the heavy scalar 
triplet Ma and the dimension-full trilinear coupling /i between the scalar triplet and 
two Higgs doublets. In this case m u ~ Y^\iv 2 jM\, where Fa is the Yukawa coupling, 
but £>r(/i — > ej) ~ Y"^M^/M A . If the scale /i is sufficiently small to suppress neutrino 
masses, Ya/Ma can be large enough to generate visible effects in rare lepton decays. A 
similar pattern can be realized also in the type-Ill seesaw, if besides a high scale Ms, a 
low scale /i, responsible for lepton number violation, is present. This has indeed been 
studied in the context of type-I seesaw [13], [10], but it can be applied here as well. In 
this case the e efl term in Eqs. fT29|) - fl3T]) is enhanced to much larger values and the x Vi 
term can be neglected. 

With such a pattern the /i — * branching ratio could be as large as ~ 10~ 4 for the 
extreme case where the Yukawa couplings would be as large as unity with triplets as 
light as few hundreds GeV. This shows that the present experimental bound is already 
relevant to exclude too large values of the Yukawas associated to too small values of 
the triplet mass. The present experimental bounds on the branching ratios give the 
following constraints on the e a p coefficients^ 

v 2 11 

| e I = lyt yi < 11 . 10 -4 / 32 ) 

: — \Y} \ \ YsU <1.5-1(T 2 (33) 
: — \Yl\^-Y^\ re <2.4-10" 2 . (34) 



I ^A*t I 



Comparison of I — > I'j and I — > 3/' decays 

The bounds of Eqs. fl32|) - fl34|) from I — > I'j decays turn out to be on the same 
parameters e as the ones obtained from \i — > 3e or r — >• 3/ decays, derived in Ref. [TO] . 
This can be understood from the fact that, at order 1/M£, for example for \i — > e7 
and /j, —>■ 3e, there is only one way to combine two Yukawa couplings and two inverse 
Ms mass matrices to induce a \i-e transition along a same fermionic line: through the 
combination e eM (i.e. the flavour structure of the /z-to-e fermionic line is the same for 
both processes, it corresponds to a \i which mixes with a fermion triplet which mixes 
with an electron). This can also be understood from the related fact that the number 
of independent parameters contained in the coefficients of the dimension five operators 
(proportional to the neutrino mass matrix) and dimension six operators (encoded in 
the e a f3 [10]) of the low energy theory (obtained in the limit of large fermion triplet 
mass) equals the number of independent parameters of the original theory. This implies 
that any physical transition studied at order 1/M|, necessarily has to be proportional 

4 Note that these bounds show that the approximation we made in the above to work only at first 
order in Y 2 v 2 /M| is justified. 
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= 1.3 


10~ 3 


• Br(n - 
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Br(r - 


■> /X7) 


= 1.3 


1(T 3 


■ Br(j - 




Br(r - 


-> e7) 


= 1.3 • 


10 -3 


■ Br(r - 


-> eee) 



to the dimension six operator coefficients, and there is only one which gives a /i-to-e 
transition: e efJi . 

As a result we obtain the following fixed ratios for these branching ratios: 

(35) 

2.1 ■ 1(T 3 • Br(r~ -> e-e + /i~) , (36) 
2.1 ■ 10~ 3 ■ Br{j- -> /i"/i + e") . (37) 

The ratios are much smaller than unity because / — > 3/' is induced at tree level through 
mixing of the charged leptons with the charged components of the fermion triplets 
|10j . while I —>■ l /r y is a one-loop process. The results of Eqs. fl35|) - fl3T|) hold in the limit 
where Ms 3> Mw,z,h, as they are based on Eq. fl3~TT) . Not taking this limit, i.e. using 
Eq. (|68|) of the Appendix, for values of as low as ~ 100 GeV, these ratios can vary 
around these values by up to one order of magnitude. Numerically it turns out that 
the bounds in Eqs. (13"2"j) - ([3~4"|) are thus not as good as the ones coming from \x —>■ eee, 
r — > eee and r — > fi/ifj, decays, which give |e e(U | < 1.1 • 10~ 6 , je^l < 4.9 • 10~ 4 , \e eT \ < 
5.1 • 10~ 4 respectively (using the experimental bounds: Br(fi — > eee) < 1 • 10~ 12 [I], 
Br{r -> eee) < 3.6 • 10~ 8 [H] and Br{r -> fififj) < 3.2 • 10~ 8 [H])Jf| This shows that 
even if the upper limits on [i — > e^ and r — > Z7 are improved in the future by three or 
two orders of magnitude respectively, the /1 — > 3e and r — >• 3/ will still provide the most 
competitive bounds on the e a p (a 7^ (3). This can be clearly seen from the bounds, 
Br(n — > ej) < 10" 15 , Br(r — > ^7) < 4 • 10~ n and _Br(r — > ej) < 5 ■ 10~ n , that one 
obtains from Eqs. (l35l) - fl3~Tj) using the experimental bounds on the I — > 31' decays. 

This leads to the conclusion that the observation of one leptonic radiative decay 
by upcoming experiments would basically rule out the seesaw mechanism with only 
triplets of fermions, i.e. with no other source of lepton flavour changing new physics. 
To our knowledge this is a unique result. 

This is different from other seesaw models. For instance, in type I seesaw, for the 
same reasons as for the type-Ill model, the ratios of Eqs. (I3~5]) - (I3~?1) are also fixed at 
order 1/Mjy, but unlike for this type-Ill model, both processes are instead realized 
at one-loop. As a result, generically, I —> l'^ dominates over I —> 31' because the 
latter suffers an extra a suppression. On the other hand, in type II seesaw, no definite 
predictions for these ratios can be done, because both types of decays depend on 
different combinations of the parameters [10] . This stems from the fact that in the 
type-II model the Yukawa coupling couples a scalar triplet to two light fermions, 
so it carries two light lepton flavour indices, instead of one in the type-I and type-Ill 
models. As a result there are several combinations of the Yukawa couplings which can 
lead to a /x-to-e transition in this mo 



5 Note that these bounds from r decays are better than the ones quoted in Tabic. 8 of Ref. [TU], as 
we have used the new experimental limits on r — > 3/ decays of Rcf . p~4j . This also leads to the new 
following bounds: |e MT | < 5.6 • 10~ 4 (from Br(r -> e+e - ^ - ) < 2.7 • 1(T 8 ) and |e er | < 7.2 • 10 -4 (from 
Br(r — > /!+ ~e~) < 4.1 • 10~ 8 ). We thank M. Nemevsek for pointing to us the existence of Ref. [14] . 

6 For instance the /.t — * 3e transition involves the combination Y&^eY^ while the /1 — * cy involve 
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5 /x to e conversion in atomic nuclei 



Beside / — > f 7 and I — > 31 decays, fermion triplets can also induce \i to e conversion 
in atomic nuclei. The relevant diagram turns out to be a tree level one, as for I — > 31 
decays, where /1 goes to e + Z with the Z connected to a u or d quark fermion line. For 
the reasons given above, or simply from the fact that this diagram involves exactly the 
same [i-e-Z vertex as the /1 — > eee decay, fi to e conversion gives a constraint on the 
same e efl parameter than from /x — > eee decay (or than from /i — > e'j decay). Using the 
experimental upper bound for the /i to e conversion rate to total nucleon muon capture 
rate ratio for ^l^i nuclei, R^ e < 4.3 • 1CT 12 [15], the bound one obtains actually turns 
out to be even more stringent than from fi —>■ eee: 

\e efl \ < 1.7 x 10~ 7 (38) 

This bound can be straightforwardly obtained by determining the quark-lepton effective 
interaction induced by the Z exchange 

c eff = -V2G F (I n a PL9^i 3 )y< 

x (u la [(l - jjsin 2 ^) - 76 )] u + d Ta [(_l + ^sin 2 ^) +75]^ (39) 
which using standard formula, for example Eq. (2.16) of Ref. [16], gives 

R ^e = l A . 1Q 1 . | £ ^|2_ ^ 4 q) 

This leads to the following fixed ratio predictions for ^Ti 

Br(fi eee) = 2.4 • KT 1 ^ 6 (41) 
Br(/2 -> erf) = 3.1 • 10^R^ e (42) 

which allows further possibilities to test and/or exclude the model. Results from the 
gold nuclei, which experimentally gives R^ e < 7 ■ 10~ 13 [17] , are of same order of 
magnitude. Note that the PRISM collaboration [18] is expected to improve the exper- 
imental bound on R^ e for the ||Tz nuclei by several orders of magnitude in the long 
term. 



6 Summary 

We have calculated the \i — > and r — > Z7 decay rates in presence of one or more 
triplets of fermions. As with right-handed neutrinos, the obtained rate is in general 
extremely suppressed but in special cases (not necessarily tuned) it can exceed the 
present experimental bounds. Unlike for other seesaw models, the observation of a 

the combination Y/^^iY^ le with I = e,[i,T see e.g. [10] . 
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leptonic radiative decay rate close to the present bounds, would nevertheless be in- 
compatible with bounds which arise in this model from / — > 31' decays. Similarly it 
would be incompatible with the bound from \x to e conversion we have determined. 
This provides an interesting possibility to exclude this model as the unique low energy 
source of lepton flavour changing new physics. 
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Appendix 



The fourteen diagrams of Fig. 1 can be grouped according to the fermion circulating 
in the loop. Performing the calculation in the 't Hooft-Feynman gauge, after loop 
integration, the various amplitudes, at 0((j^) 2 ), are: 



rp<t>-,W- 



II - 



+ T, 



.G 



SM 



; m^u e (p - q) (1 + 7s) (2p • e) u„ (p) 



y/2 32tt 2 

+ (eU 0l/ „) ei (ul) F 2 (x Vi ) + (U 0u „] 
\ /in 



I ft 



(43) 



,W~ rpcj>- _|_ rp<t>-,W- _|_ rpW 



-w- 



.G S F M e 
' V2 32tt 



2 m^T e (p - q) (1 + 7s) (2p ■ e) u„ (p) 



{ (^^ei + { Y ^) et ( M ^ 6 )i, J^M^ 



M 2 wl 



M 2 W 



2 \ 



(44) 
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Z,H,rj 



Tl + T«+Tl 
r<SM 

Mr? 



-m^Ue (p - q) (1 + 75) (2p ■ e) u M (p) 

,2 



v^2 327T 2 

2 



IT 



[ F^M^ 1 ) . (M^Y^e) . (F 9 fa) + F 9 (z Si )) 



(45) 



T, 



Z,H,r] 



Tf + T, H + TP 



r<SM 

. <-J J? 



V2 32tt 2 

_ m % 

G(x) are the following functions: 



"fy«e (p -?)(! + 7s) (2p • e) m m (p) e eM G (j/h, , (46) 



where x Vi = x Ei = ^ = ^, 



mf m 



^i(x) 

^ 3 (x) 
F 4 (x) 

F 7 (x) 
F s (x) 
F 9 (x) 



10 - 43a; + 78x 2 - 49x 3 + 4a; 4 + 18a; 3 log(x) 
3(-l +x) 4 

2(5 - 24x + 39a; 2 - 20x 3 + 6x 2 (-l + 2x) log(x)) 
3(-l +x) 4 

7 - 33x + 57a; 2 - 31x 3 + 6x 2 (-l + 3x) log(x) 
3(-l +x) 4 

-38 + 185a; - 246a; 2 + 107a; 3 - 8a; 4 + 18(4 - 3x)x 2 log(x) 

3(-l + x) 4 
1 — Qx + 3a; 2 + 2a; 3 — 6x 2 log(x) 

3(-l +x) 4 
7 - 12x - 3x 2 + 8x 3 - 6x(-2 + 3a;) log(x) 
3(-l +x) 4 

40 - 46a; - 3a; 2 + 2x 3 + 7x 4 + 18x(4 - 3x) log(x) 
3(-l +x) 4 

x(-16 + 45a; - 36a; 2 + 7a; 3 + 6(-2 + 3x) log(x)) 

3(-l +x) 4 
x(2 + 3a; — 6a; 2 + x 3 + 6a;log(x)) 
3(-l +x) A 



and Fi(x) and 

(47) 
(48) 
(49) 
(50) 
(51) 
(52) 
(53) 
(54) 
(55) 
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+ 5^ 



- - cos 2 9 



4 - 9y k + 5yl + 6(1 - 2y h )y h log(^ ; 



2/, 



6(-l+Z//J 4 
16 - 45^ + 36^ - H ~ 6 (- 2 + 3z h) log(^) 



+8 ( ^ - cos 2 6W 



-8 (1 - cos 2 w ) 



2(-l + ^) 4 
4 - 91/,, + 5yf a + 6(1 - 2y h )y h \og{ yi 

6(-l + ^) 4 
2(-l + y 2 -2^1og(^)) 



-yii 



-20 + 39y h - 24y 2 + 5^ + 6(-2 + y,J log(y h ] 
6(-l +yzj 4 



(56) 



Since y\ v Zi v x Vi <C 1, it is a good approximation to take the lepton flavour conserving 
quantities y\ i and Zi { to zero and to keep only the linear term in the flavour changing 
quantities x n : 



(57) 
(58) 





10 

3 






10 

3 


8 


F 3 (x Ui ) 


7 
3 


5 

3^ 




= C = 


-6,56 



(59) 
(60) 



Summing over i and neglecting terms of 0((Y^v /Ms) n ) with n > 2, we obtain: 



Tf 



i 

.Or p 



V2 32tt 2 
7 



m^u e (p-q)(l+ 75) (2p • £) m m (p) 



/; 



PMNS) 



U 



PMNS / . 



(61) 



,w- 



.Gp M e 



;m^u e (p - g) (1 + 75) (2p ■ e) (p) 



327T 2 



(62) 



' V^2 32vr 2 



iTiu'U'e (p - 9) (1 + 75) (2p • e) u„ (p) e eA1 x C (63) 
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rpZ.H,r] rpZ,H,T) 



.Gp M e 



i —r^^r^ m ^ u e (P -?)(! + 75) (2p ■ e) Uf, (p) 



y/2 32tt 2 

{y e ^ + E y (^ M ^L ( M s ly =)* + } 



(64) 



where 



x 



-30 + 153x E . - 198x^ +75xf, +18(4-3x s .)a;^ logx S - , 

A(xsJ = k * * y SJ s ' s ' (65) 

u( , 33 - 1% S< - 45y| + 30^ + 18(4 - 3y s< )y Si logy Sl 

= (66) 

C{Z ^ = H^W • (67) 

The total amplitude is then: 

qsm e 

T(ji^>ery) = i F 2 m M w; (p - q) (1 + 75) (2p ■ e) ?v (p) (68) 
( y + C J e eM - ^ x Vi {UpMNs) ei (u PMNS ^ _ + 

E y (^ M ^) e , ( M £ ly ^ + + ^ 

i 

This result is valid at 0((^^) 2 ). For x^^y-E^ z^ ^> 1, the additional limit x^, 2/^, z^. 
00 can be taken, which leads to the result displayed in the text, Eq. 029|) . 
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